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NAVIGATION FOR lUS DEPLOYMENT 

TDRSS NAVIGATION ACCURACY IN SUPPORT OF lUS DEPLOYMENT PHASE 1 

By Alan D, Wylie 
Mathematical Physics Branch 


SUMMARY AND INTRODUCTION 


The first phase of a study has been completed to determine the navigation accu- 
racy for tracking the Orbiter prior to interim upper stage (lUS) deployment 
using the tracking data relay satellite system (TDRSS), The study was performed 
to examine the Orbiter navigation accuracy for both one and two TDRSS satellites, 
for short and long data arcs, and for Doppler-only and Doppler-plus range sol- 
utions. All test cases were run with the Orbiter in a 150-n. mi. circular 
orbit, 28.5 degree inclination, at the time interval from the completion 
of the orbital maneuvering system (OMS)-2 maneuver to OMS-2 plus 2 hours 
(approximate time for lUS deployment). The data used were simulated by the 
simulation navigation (SIMN.AV) program (ref. 1), The software tool used 
to process the TDRS data was the Shuttle Navigation Analysis Program (SNAP) 
(ref.l), a Kalman filter tool used to solve for the Orbiter position and 
velocity. Results of the study were used in a presentation to NASA Head- 
quarters to summarize the expected navigation accuracy using the TDRS syitem. 

The study is now being continued to analyze in more detail the effects of 
length of data arc, type of data, and number of TDRS satellites and to exa- 
mine the effects of vehicle venting, measurement noise, and TDRS ephemeris 
ert ors on navigation solutions in support of lUS deployment. 

Conclusions drawn from the first phase of the TDRS navigation study are as 
follows . 

a. Data from both TDRS satellites were essential for accurate navigation 
results . 

b. Range data were essential for the short arc test case (At =8 min) but 
were not needed for the long arc test case (At = 55 min). 

c. With Doppler and range data from both TDRS satellites, the results con- 
verged to a reasonable solution after 5 to 10 minutes of data. 

Further analysis in these three areas and in the effects of vehicle venting, 
measurement noise, and TDRS ephemeris errors is now in progress. 
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DE SCIUPTION OF TEST CASES 

Test cases were defined to detennine the effects on navigation accuracy by the 
fol lowing; 

a. Number of TORS satellites 

b. Doppler only versus Doppler plus range solutions 

c. Length of tracking arc 

Table I summarizes the eight test cases that were run. Four of the cases (lA 
to 4A) used both Doppler and range data as input; the other four cases (IB to 
40J processed only Doppler data. Cases lA and IB are the long data arc cases 
with data from both satellites. Cases 2A and 2B are the long arc, one TORS 
cases. Cases 3A and 3B are the short arc cases. Cases 4A and 4B are the same 
as lA and IB, except that the position of the tracking arc is shifted to allow 
data from both TORS satellites to be available early in the tracking arc. 

The basic orbit used for the test cases was from the Department of Defense 
vDOD) geosynchronous mission, which is one of the six basic reference missions 
being used in IDS studies (ref. 2). The data tracking arcs were positioned be- 
tween the end of the OMS-2 burn and OMS-2 plus 2 hours, which is the approximate 
time of the reference mission for lUS deployment. In this time frame for the 
ODD mission, the Orbiter is in a 150-nautical-mlle orbit, 28.5-degree inclina- 
tion. The two TORS satellites are positioned at longitudes of 189 E and 319 E, 
which are the baseline positions now being assumed for all studies involving 
TORS. 

Fable II summarizes the inputs for the eight test cases. Parts A through E 
contain the inputs common to all runs. Part F includes the parameters that 
varied among the eight runs. 


DESCRIPTION OF ERROR MODEL 

Errors applied in the study included bias and noise to the Doppler and range 
data, TORS ephemeris errors, and a continuous 0.5-pound vent along the vehi- 
cle's velocity vector. Table Ilb summarizes the actual values used. 


FILTER WEIGHTING INPUTS 

Weighting factors for SNAP program input were set to inform the filter of the 
known perturbations being applied in the simulation of the input data. These 
’weighting factors include the program's uncertainty for data noise and bias, 
and a state noise constant (SNC) modeled into the program to allow for Orbital 
perturbations known to exist but not modeled into the program (such as drag 
perturbations, venting, etc.). The SNC value was set to allow for the propa- 
gation error due to the venting. (A 0.5 pound vent contributes to a 0.5 nauti- 
cal mile downtrack propagation error after one revolution.) The filter uncer- 
tainties for Doppler noise and bias and for range noise were set to the values 
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of ih(3 noise and bias being applied to the data. The filter uncertainty for 
range bias was set to allow for the range bias that was applied to the data 
and for the bias that was introduced by the radial position error simulated 
for the two TORS satellites. The filter weighting values are summarized in 
table lie. 


RESULTS AND CONCLUSIONS 


In this study the procedure for evaluating the navigation accuracy for each 
run was to set the filter weights (estimates of measurement noise and bias 
and the state noise constant) to account for the known errors in the data 
(measurement noise and bias, TDRS ephemeris errors, vehicle venting) and ex- 
amine the filter output of its uncertainty in Orbiter position and velocity 
in the local-vertical (radial, downtrack, crosstrack) coordinate system. 

Tables III and IV contain tabulations of the filter uncertainty for each com- 
ponent of Orbiter position and velocity for critical times in each of the 
eight test cases. Figures 1(a) through 8(b) show time history plots of these 
values for the eight runs. Tables V through VII provide interesting compari- 
sons which have led to the following conclusions. 

a. Data from both TDRS satellites were essential for accurate na'igation 
results. 

b. Range data were essential for the short arc test case (At = 8 min) but 
were not needed for the long arc test case (At = 55 min). 

G, Short arc solutions were acceptable only when processing Doppler and range 
data from both TDRS satellites. 

Another interesting comparison is between the results for the best geometry 
case (case lA), with two TDRS satellites and Doppler and range data, against 
the currently documented TDRS navigation accuracy numbers obtained from refer- 
ence 3. This comparison is shown in table VIII as a point of interest. Note 
that the in-plane results from case lA are not as accurate as the currently 
documented numbers and that the out-of-plane results are an improvement to the 
prior results. However, case lA results are based on simulated data and derived 
from an error model which is being further analyzed and should therefore be 
regarded as preliminary. The future results will be documented in a report to 
be published. 

Further work is necessary prior to determining the final navigation accuracy 
numbers in support of lUS initialization. Test cases with orbital inclinations 
of 37.4 degree and 57 degree will be analyzed. Monte Carlo studies will be 
performed using simulated data and different error models. Finally, a review 
of real data results from the Apollo/Soyuz Mission will be made in an attempt 
to evaluate the error model being used for the simulated data results. 
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TABLE I.- DESCPJPTION OF TEST CASES 


Case 

number 

Type of 
data 

Number of 
TORS 

satellites 

Length of 
tracking arc, 
min 

Length of 
propagation 
interval , 
min 

lAa 

Doppler + range 

2 

55 

60 

Ik 

Doppler + range 

1 

49 

60 

3A 

Doppler + range 

2 

7.7 

60 

4Ab 

Doppler + range 

2 

67 

0 

1B3 

Doppler only 

2 

55 

60 

2B 

Doppler only 

1 

49 

60 

3B 

Doppler only 

2 

1.1 

60 

4Bb 

Doppler only 

2 

67 

0 


aSecond TORS not available until 40 min elapsed in tracking arc, 
^Both TDRS's available within first 4 min of tracking arc. 
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TABLE II.- PROGRAM INPUTS 

a. State vectors, RNP matrix, Greenwich hour angle 

(1) Time tag for Orbiter and TORS state vectors 
G.m. t. - March 30, 1982 13 hr 44 min 45.574 sec 
Get - 0 hr 44 min 45.574 sec 

(2) Orbiter state vector (in and m/sec) 

X -1242.2912640 + 03 X -7278.74653 

Y 5838.248874D +03 Y -2246.207141 

Z -2945.056326D +03 Z -1382.356544 

(3) TORS State vector - longitude “ 1890 (m m/sec) 

X -31134159.96 X 2073.412485 

Y -28433920.33 Y -2270.336867 

Z 95830.30872 Z -6.4870 

(4) TORS state vector - longitude « 319<> (m and m/sec) 

X 41794192.34 % 406.40685 

Y -5573291.5 Y 3047.67829 

Z -129746.227 Z -1.18327 

(5) RNP matrix 

-.993267275 -.115804370 .00307723788 

.115803735 -.993272041 -.000384193237 

.00310102560 -.0000252508913 .999995191 

(6) Greenwich hour angle - 187.05880 
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TABLE II.“ Continued 


b. Errors modeled Into study 
U) TORS epheinerls en ors 

Position; radial 15 m® 
downtrack 66 mb 
crosstrack 6' m 

Velocity: radial .005 m/secb 

downtrack .0005 m/sec 
crosstrack 0 

(2) Measurement errors 


Doppler noise .004 Hz 

Doppler bias 0 

Range noise 3 in 

Range bias 3 in 


{3} Orbiter venting 


Magnitude 0.5 pounds continuous along vehicle velocity vector 

Orbiter cross-sectional area perpendicular to velocity vector 
249.909 m2 


Orbiter mass 94345.30707kg 
c. Filter weights 

(1) Variance for Doppler bias 

(2) Variance for Doppler noise 

(3) Variance for range bias 

(4) Variance for range noise 

(5) State noise constant 

(6) All time constants 


0 

0.16D-4 Hz 2 
909 m2 
9 m^ 

1 .603D-06 m2/sec3 
400 sec 


^The Earth's gravitational constant, pg was changed accordingly in computing 
the gravitational effects on the TORS satellites. 


^Error introduced into the two TORS satellites by clianging the time tag of the 
reference state vectors by 0.021 sec. 
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TABLE II.- Continued 


(7| Weighting factor control non! Inearltles 1.2 

in the measurement incorporation equa- 
tions 

(B) Minimum elevation above Earth horizon 0.2o {/g.s n* mi.) 

for processing TDRS measurements 

d, Initial errors in Orbiter state vector 


(1) 

Position: 

radial 

downtrack 

crosstrack 

45/1.65 m {15 000 ft) 

14 635.57 m (48 000 ft] 
3354.10 m (11 000 ft) 


(2) 

Veloci ty: 

radial 

downtrack 

crosstrack 

16,16 m/sec 
3.96 m/sec 
3,35 m/SBC 

(53 ft/sec) 
(13 ft/sec) 
(11 ft/sec) 


Initial filter 

state covariance matrix 



(1) 

Posi tion: 

radial 

downtrack 

crosstrack 

3 X 104 ,n 
3 X 10^ m 
3 X 104 m 



(2) 

Velocity: 

radial 
feiii^rack 
t track 

31,6 m/sec 
31,6 m/sec 
31,6 m/sec 



Start and stop 

times 




(1) 

Filter 

processing 

Cases 
lA, IB 

Cases 
2A, 2B 

Cases 
3A, 3B 

Cases 
4A, 4B 


Start 

13:44:45,5/4 13:44:45.5/4 13:44:45.574 

13:44:45.5/4 


Stop 

15:44:45.! 

574 15:44:45.574 15:44:45.574 

15:44:45.574 


Delta, min 

120 

120 

120 

120 
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TABLE II.- Concluded 


Data 

Inter^/al 

Cases 

lA, IB 

Cases 
2A, 2B 

Cases 
3A, 38 

Cases 
4A, 4B 

Start 

13:49:52.7 

14:10:50 

14:36:48 

14:35:48 

Stop 

14:44:45.574 

14:44:45.574 

14:44.45.574 

15:44:45.574 

Delta, min 

55 

34 

7.7 

67.7 

Orbiter 

venting 

interval 

Cases 
lA, 18 

Cases 
2A, 2B 

Cases 
3A, 3B 

Cases 
4A, 4B 

Start 

13:44:45.574 

13:44:45.574 

13:44:45.574 

13:44:45.574 

Stop 

15:44:45.574 

15:44:45.574 

15:44:45.574 

15:44:45.574 

Delta, tnin 

120 

120 

120 

120 
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TABLE m.- TABULATION OF FILTER COVARIANCE OUTPUT 


AT TIMES OF INTEREST; INPUT DATA - DOPPLER AND RANGE 


Time, 

sec 

Comment 

ft 

°V. 

ft 

ft 


-V, 

fps 

-W' 

fps 

Run la 

- Two TDRS satellites; 

best geometry 





5HOH5 

Prior to acquisition 
of second TDRS 

286 

2976 

16 074 

3.3 

0.14 

13.47 

52055 

After acquisition 
of second TDRS 

90 

63 

271 

0.16 

0.08 

0.17 

53085 

End of data arc 

66 

57 

104 

0.09 

0.09 

0.15 

56685 

End of propagation 

648 

1866 

186 

2.0/ 

0.55 

0.23 

Run 2a 

- One TDRS satellite only; poor geometry 





52815 

End of data arc 

377 

1449 

8 968 

1.45 

0.20 

3.01 

56685 

End of propagation 

729 

5750 

832 

6.28 

0.68 

10.83 

Run 3a 

- Two TDRS satellites; 

short data 

i arc 





53085 

End of data 

97 

107 

238 

0.54 

0.49 

0.81 

56685 

End of propagation 

1591 

4949 

662 

5.22 

1.39 

0.44 

Run 4a 

- Two TORS satelV'es; 

data available from both 

satellites 

early in 



tracking arc 







52815 

TDRS 1 LOS 

(after 3 min of data) 

168 

95 

305 

1.9 

0.75 

2.2 

549^5 

TDRS 2 LOS 

220 

86 

264 

0.37 

0.07 

0.19 

55615 

TDRS 1 AOS 

430 

403 

174 

0.79 

0.32 

0.26 

56685 

End of tracking arc 

75 

125 

172 

0.10 

0.13 

0.20 
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TABLE IV.- TABULATION OF FILTER COVARIANCE OUTPUT 



AT TIMES OF INTEREST; 

INPUT DATA 

- DOPPLER 

ONLY 



Time , 
aec 

Comment 

<7u. 

ft 

C-v» 

ft 

Oyt 

ft 

fps 

fps 

Jw» 

fps 

Run 1b 

- Two TORS satellites 

beat geometry 




52025 

Prior to TORS 2 AOS 

ll6i| 

5 612 

34 965 

5.69 

0. 18 

21 .93 

52055 

After TDRS 2 AOS 

m 

1 138 

7 995 

0.98 

0.19 

3.30 

53085 

End of data 

102 

8i| 

168 

b.12 

0.11 

0.24 

56685 

End of propagation 

650 

1 936 

200 

2.15 

0.55 

0.32 

Run 2b 

- One TDRS satellite! 

poor geometry 





32815 

End of data 

2 086 

18 06i| 

93 153 

20.85 

1 .06 

19.91 

56685 

End of propagation 

2 962 

37 401 

2 666 

42.87 

1.98 

110.25 

Run 3b 

- Two TDRS satellites! 

short data arc 





1 53085 

End of data 

1 972 

1 743 

5 440 

2.36 

2.41 

3.24 

56685 

End of propagation 

11 332 

58 234 

4 907 

61.83 

8.53 

4,17 


Run 4b - Two TORS satellites; d ata available from both satellite? early in 
bracking arc 


52815 

TDRS 1 LOS (after 
3 rain of data) 

4 741 

2 432 

13 094 

8.90 

3.42 

11.54 

54965 

TDRS 2 LOS 

294 

413 

946 

0.72 

0.15 

0.55 

55615 

TDRS 1 AOS 

109 

355 

867 

0,35 

0.06 

0.47 

56685 

End of tracking arc 

117 

129 

198 

0.13 

0.14 

0.29 


OF l?OOR 
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TABU V.- TWO TORS VERSUS ONE TORS SOLUTIONS 
(LONG TRACKING ARC) 


Component 

Doppler 

plus range 

Doppler 

only 


Two 

TORS 

results 

One 

TORS 

results 

Two 

TORS 

results 

One 

TORS 

results 

Oyi , f t 

648 

729 

650 

2962 

Ov* 

1866 

5750 

1936 

37 401 

, ft 

186 

3832 

200 

2666 

0 ., fps 
u 

2.07 

6.28 

2.15 

42.87 

0 ., fps 

V 

0.55 

0.68 

0.55 

1.98 

0 ., fps 
w 

0.23 

10.83 

0.32 

110.25 


^Not stable. 


Conclusions; 

(1) For Doppler plus range data, two TDRS’s are necessary to get good plane 
estimate. 

(2) F-ir Doppler only, two TDRS's are necessary for both in-plane and out- 
of-plane estimates. 


12 



TABU VI.- DOPPLER PLUS RANGE VERSUS DOPPLER ONLY {TWO TORS' S) SOLUTIONS 


Component 

Long 

(At 

tracking arc 
= 55 min) 

Short tracking arc 
(At = 8 min) 


Doppler 

only 

Doppler 

plus 

range 

Doppler 

only 

Doppler 

plus 

range 

ou> 

650 

648 

11 332 

1591 

Oy > ft 

1936 

1866 

58 234 

4949 

O' W > ^ 

200 

186 

4907 

662 

a,, fps 
u 

2.15 

2.07 

61.83 

5.22 


0.55 

0.55 

8.53 

1.39 

0 > j fps 
w 

0.32 

0.23 

4.17 

0.44 


Conclusions : 

(1) Doppler only sufficient for the long tracking arc. 

(2) Range data essential for the short tracking arc. 
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TABLE VII.- LONG ARC VERSUS SHORT ARC SOLUTIONS (TWO TORS'j) 


Component 

Doppler plus 

range 

Doppler only 


Long arc 
resul ts, 
(At ” 55 min) 

Short arc 
results, 
(At = 8 min) 

Long arc 
results, 
(At = 55 min) 

Short arc 
results 
(At = 8 min) 

Oj I ft 

648 

1591 

650 

11 332 

ay, ft 

1866 

4949 

1936 

58 234 

0\f) f ft 

186 

662 

200 

4907 

0 .. fps 
u 

2.07 

5.22 

2.15 

61.83 

0 ., fps 

V 

0.55 

1.39 

0.55 

8.53 

0 ., fps 
w 

0.23 

0.44 

0.32 

4.17 

Concluding 

remark: 





Short ?rc solution sufficient only when processing Doppler plus range from 
two TMSS satellites. 
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TABLE VIII.- CASE U RESULTS VERSUS 


CURRENT OFFICIAL TDRS NAVIGATION ACCURACIES (3a) 


Component 

Current TDRS 
aeouraoy numbers 

Case sA resultL. 

u, ft 

1000 

1950 

V, ft 

3000 

5600 

w, ft 

3000 

560 

u, fps 

3.0 

6.0 

V, fps 

1.0 

1 .b 

w, fps 

3.0 

.7 


\ 


Position error 



(%) Position. 

Figure 1.- Navigation accuracy with TDRS, case 1A. 
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Velocity error 1°, fps 



(b) Velocity. 
Figure 1.- Concluded. 
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Position error lo, ft 



0 20 UO 60 80 100 120 

Time from 0KS>2 cutoff, min 


(a) Position. 

Figure 2.- Navigation accuracy with TORS, case 2k. 
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^ Downtrack 
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(b) Velocity. 
Figure 2.- Concluded. 
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(a) Poaltlon. 

Figure 3 .- Navigation accuracy with TORS, case 3A. 
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Figure 3.- Concluded. 
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Position 



Time from OMS-2 outoff, min 


(a) Position. 

Figure 4.- Navigation accuracy with TDRS, cise 4A. 
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(b) Velocity. 
Figure Concluded. 
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Figure 5.- Navigation accuracy with TORS, case IB 
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(b) Velocity. 
Figure 5.- Concluded. 
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Figure 6.- Navigation accuracy with TORS, case 2B. 
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Figure 7.- Navigation accuracy with TORS, case 3B. 
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Figure 8.- Navigation accuracy with TORS, case UB. 
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Figure 8.- Concluded. 
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